#######第一题
m<-m[,-1]
m<-m[,-1]
m<-m[,-33]
m<-as.matrix(m)
hist(log2(m))
#######第三题
install.packages("ggplot2")
library("ggplot2")
up <- prostat[prostat$FC > log2(1.2) & prostat$P < 0.05, ]
do <- prostat[prostat$FC < -log2(1.2) & prostat$P < 0.05, ]
ggplot(data = prostat, mapping = aes(x = FC, y = -log10(P))) +
  geom_point(color = "black", alpha = 0.5) +  # 所有数据点
  geom_point(data = up, mapping=aes(color = "red")) +  # 上调蛋白
  geom_point(data = do, mapping=aes(color = "blue")) +  # 下调蛋白
  labs(x = "log2 Fold Change", y = "-log10(P)", title = "Volcano Plot")
#######第二题
m <- t(m)

# 对转置后的矩阵进行层次聚类分析
clustering_result <- hclust(dist(m))

# 画出树状图
plot(clustering_result)
